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Chemical bonding at the active site of hen egg-white lysozyme 
(HEWL) is analyzed on the basis of Bader's quantum theory 
of atoms in molecules [QTAIM; Bader (1994), Atoms in 
Molecules: A Quantum Theory. Oxford University Press] 
applied to electron-density maps derived from a multipole 
model. The observation is made that the atomic displacement 
parameters (ADPs) of HEWL at a temperature of 100 K are 
larger than ADPs in crystals of small biological molecules at 
298 K. This feature shows that the ADPs in the cold crystals of 
HEWL reflect frozen-in disorder rather than thermal vibra- 
tions of the atoms. Directly generalizing the results of 
multipole studies on small-molecule crystals, the important 
consequence for electron-density analysis of protein crystals is 
that multipole parameters cannot be independently varied in a 
meaningful way in structure refinements. Instead, a multipole 
model for HEWL has been developed by refinement of atomic 
coordinates and ADPs against the X-ray diffraction data of 
Wang and coworkers [Wang et al. (2007), Acta Cryst. D63, 
1254-1268], while multipole parameters were fixed to the 
values for transferable multipole parameters from the 
ELMAM2 database [Domagala et al. (2012), Acta Cryst. 
A68, 337-351] . Static and dynamic electron densities based 
on this multipole model are presented. Analysis of their 
topological properties according to the QTAIM shows that the 
covalent bonds possess similar properties to the covalent 
bonds of small molecules. Hydrogen bonds of intermediate 
strength are identified for the Glu35 and Asp52 residues, 
which are considered to be essential parts of the active site of 
HEWL. Furthermore, a series of weak C— H- ■ O hydrogen 
bonds are identified by means of the existence of bond critical 
points (BCPs) in the multipole electron density. It is proposed 
that these weak interactions might be important for defining 
the tertiary structure and activity of HEWL. The deproto- 
nated state of Glu35 prevents a distinction between the 
Phillips and Koshland mechanisms. 
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1. Introduction 

The three-dimensional structure of the antibacterial hen egg- 
white lysozyme (HEWL) was uncovered in 1965 (Blake et al, 
1965), heralding enzyme crystallography. Today, lysozyme 
resides among the most intensively studied enzymes and its 
function is well known. Lysozyme damages bacterial cells by 
cleaving the yS(1^4) glycosidic hnkage between alternating 
units of N-acetylmuramic acid and A/^-acetylglucosamine, 
which are the building blocks of bacterial cell walls. According 
to the Phillips mechanism (Philhps, 1966), the two residues 
Glu35 (glutamic acid) and Asp52 (aspartic acid) (Fig. 1) 
play an essential role. The terminal proton of Glu35 is trans- 
ferred to the O atom of the glycosidic bond between two 
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neighbouring sugar residues, leading to the cleavage of the 
glycosidic bond and the formation of a carbenium ion. The 
positive charge of this carbenium ion is stabilized by the 
negative charge of Asp52 until a hydroxyl ion binds to the 
positive C atom and Glu35 is reprotonated. 

Another reaction mechanism, proposed by Koshland 
(1953), involves a glycosyl-enzyme intermediate in which the 
substrate is covalently bonded to Asp52. Vocadlo et al. (2001) 
concluded in their electrospray-ionization mass-spectrometric 
study that the catalytic mechanism of lysozyme involves an 
intermediate complex in which one sugar ring is covalently 
bonded to Asp52. A crystallographic study at 1.5 A resolution 
of HEWL crystallized with substrate by Strynadka & James 
(1991) supports the Phillips mechanism (Philhps, 1966) and 
does not indicate a covalent intermediate. They pointed out 
the environment of Glu35 and Asp52, which among other 
things is constructed by hydrogen bonds that can assist the 
catalysis according to Phillips (1966). 

An important factor for the catalytic or other functions of 
a protein is its flexibility (Branden & Tooze, 1999). All atoms 
are moving in a temperature-dependent breathing motion, 
which is retained in crystals. In crystal structures, flexible parts 
appear as disorder or give rise to large values of B factors, 
while these regions attain low levels of electron density or 
almost no density at all (Branden & Tooze, 1999). Analysis of 
crystallographic B factors thus allows inferences concerning 
the relations between structure and dynamics of proteins 
(Frauenfelder et al, 1979; Artymiuk et al, 1979; Karplus & 
Schulz, 1985; Carugo & Argos, 1999; Parthasarathy & Murthy, 
2000; Stocker et al, 2000; Radivojac et al, 2004; Yuan et al, 
2005; Weiss, 2007). 

Apart from flexibihty, the stability and the function of 
proteins primarily depends on chemical interactions. The 
electron density in the unit cell can provide quahtative and 
quantitative information on chemical interactions of mole- 
cules, including proteins (Gatti & Macchi, 2012). This 
encompasses information on the degree of covalency of bonds, 
on the nature of hydrogen bonds and on covalent and 




Figure 1 

Schematic drawing of Glu35 and Asp52 within HEWL. 



electrostatic interactions, opening the door for a detailed 
insight into protein structure, function and potentially to drug 
design. 

Electron-density studies are an established tool for 
extracting information about chemical interactions involving 
small molecules, and thus can shed light on the functionality 
regarding chemical stability, chemical reactions and physical 
properties. Electron-density studies are usually based on 
refinement of the parameters of the multipole model against 
X-ray diffraction data (Hansen & Coppens, 1978). The 
multipole model is an extension of the independent atom 
model (lAM) by many functions and many parameters that 
describe the rearrangements of electrons as occur, for 
example, upon the formation of chemical bonds (Coppens, 
1997). Refinements of multipole parameters are only possible 
if accurate X-ray diffraction data have been measured at low 
temperature {T < 100 K) and to high resolution (<iniin> 
preferably below 0.50 A). A single protein (crambin; Jelsch 
et al, 2000; Schmidt et al, 2011) diffracts to this limit, while 
subatomic resolution data sets of several more proteins have 
been measured (Ko et al, 2003; Howard et al, 2004; Wang et 
al, 2007). Here, we present an analysis of the d^^^ = 0.65 A 
data set of HEWL as pubhshed by Wang et al (2007). 

The current and developing state of the art of X-ray 
diffraction promises more high-resolution, high-quality data 
of proteins in the future (Duke & Johnson, 2010; Dauter et al, 
2010; Mueller et al, 2012; Garman & Weik, 2013). One factor 
limiting the number of very high-resolution data sets in the 
Protein Data Bank (PDB; Berman et al, 2000) is the geometry 
of many diffractometers for protein diffraction. With an area 
detector placed perpendicular to the beam, and a maximum 
angle of incidence onto the detector of 60° for good-quality 
data, the resolution is given by the wavelength of the radia- 
tion, dmin = ^, and is usually worse than 0.50 A, even if the 
crystal could diffract better. Therefore, diffractometers with 
an area detector that can be rotated to offset positions in 26 
are highly desirable to obtain high-quahty diffraction data at 
subatomic resolutions. 

Information on chemical interactions is contained in the 
topological properties of static electron densities, where the 
latter are electron densities after the removal of all thermal 
motion. They are obtained from a structure model with all 
atomic displacement parameters (ADPs) set to zero, resulting 
in the deconvolution of electron density and atomic displa- 
cements. Topological properties of static electron densities are 
usually interpreted according to the quantum theory of atoms 
in molecules (QTAIM; Bader, 1994). This method has been 
extensively applied to crystals of amino acids and small 
peptides, especially tripeptides (Jelsch et al, 1998; Destro et al, 
2000; Flaig et al, 2002; Rodel, 2003; Scheins et al, 2004; 
Dittrich et al, 2005; Mebs et al, 2006; Checihska et al, 2006; 
Kalinowski et al, 2007; Johnas et al, 2009). As opposed to 
static electron densities, dynamic electron densities employ 
both the atomic coordinates and the atomic displacements 
resulting from a structure refinement against diffraction data. 
Thus, dynamic densities describe a smeared electron-density 
distribution. Here, we consider dynamic model densities. 
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which are computed by inverse Fourier transformation of a 
complete set of structure-factor amplitudes and phases that 
have been calculated from a structure model up to a resolution 
far beyond experimental access (Mondal et al, 2012). As a 
consequence, dynamic model densities are not subject to 
series-termination effects of the Fourier transform, thus 
allowing an accurate description of its fine topological prop- 
erties. The dynamic model densities can be accurately 
computed for any structure model, including lAM and 
multipole models. They are essentially different from the 
dynamic densities sometimes defined as the inverse Fourier 
transform of the phased experimental structure factors. 
Although the latter are closer to the experiment, the limited 
resolution of the experimental data necessarily leads to series- 
termination effects in the corresponding electron density, thus 
preventing a meaningful analysis of topological properties. 

An increasing size of the molecule is accompanied by an 
increasing number of parameters to be refined. The assign- 
ment of multipole parameters to a protein structure boosts the 
number of parameters even further, which may limit the 
refinement of the multipole model. However, it has been 
found that the multipole parameters of an atom in a particular 
environment (e.g. the C" atom in alanine) are nearly the same 
in all compounds containing this atom (Pichon-Pesme et al, 
1995). Based on this principle, several databases of transfer- 
able multipole parameters have been developed (Pichon- 
Pesme et al, 1995, 2004; Volkov et al, 2004; Dittrich et al, 2006; 
Zarychta et al, 2007; Dominiak et al, 2007, 2009; Domagala et 
al., 2012). They offer the possibility of constructing a multipole 
model of any molecule without the need to refine parameters 
beyond those of the I AM, i.e. only atomic coordinates, ADPs 
and site occupancies are refined. Electron densities based on 
multipole models of proteins have previously been considered 
for crambin (Jelsch et al, 2000; Schmidt et al., 2011), human 
aldose reductase (Guillot et al, 2008) and trypsin (Schmidt 
et al., 2003). Multipole parameters have been refined in these 
cases, whereas for scorpion toxin (Housset et al, 2000) a 
multipole model was developed with multipole parameters 
fixed to database values. A similar approach was used for a 
series of neuraminidase-inhibitor complexes (Dominiak et al, 
2009). 

The original aim of the present study was to elucidate the 
reaction mechanism catalyzed by HEWL from consideration 
of electron-density distributions. We have not achieved this 
goal yet, mainly because high-resolution X-ray diffraction 
data are only available for native HEWL (Wang et al, 2007) 
and not for protein-substrate complexes. 

However, in the course of this study we have obtained 
several results that we consider important for future electron- 
density studies of HEWL and of proteins in general, and 
potentially of other large biomolecules such as DNA or RNA. 
We present here a multipole model for HEWL. We make the 
observation that the ADPs of atoms in the ordered part of 
HEWL have values comparable to or larger than those of 
serine at 298 K, despite the fact that the diffraction data were 
measured at T = 100 K. We explain this feature as frozen-in 
disorder, but the implication is that ADPs are too large for a 



Table 1 

Crystallograpliic data of lien egg-white lysozyme (PDB entry 2vbl; Wang 
et al, 2007). 



Space group PI 

Z 1 

a (A) 27.07 

b (A) 31.25 

c (A) 33.76 

a (°) 87.98 

fi 1°) 108.00 

Y {") 112.11 

V{A^) 25057.0 

f(OOO) 10268.0 

T (K) _ 100 

Wavelength (A) 0.65 

[sin(e)/A.]„,, (A-') 0.77 

d™„ (A) 0.65 

pH 4.7 

Multiplicity 7.1 

Emerge (%) 4.00 

Completeness (%) 97.6 

Unique reOections (observed/all) 166259/187165 

«all (%) 8.39 



lAM before elongation of bond lengthsf 

Rf [fobs with |/„b!, - 4aic|/cf(/<,b») < 45/all reflections] (%) 8.39/8.54 

^wF [fobs with |/„bs - 4aic|/o-(4bs) < 45/all reflections] (%) 10.19/10.80 

Final lAMf 

-Rf [fobs with - /caic|/'j(4bs) < 45/aU reflections] (%) 8.39/8.54 

^wF [fobs with j/obs - /caic|/<^(4bs) < 45/all reflections] (%) 10.19/10.80 
ELMAM2 modelt 

i?F [fobs with \U, - /„,i,|/ff(/obs) < 45] (%) 8.07 

^wF [fobs with |/„bs - /„,ic|/ff(4bs) < 45] (%) 9.88 



t Agreement indices from the lAM and the ELMAM2 model are from the present 
work. 

meaningful refinement of the multipole parameters (Dittrich 
et al, 2009; Gatti & Macchi, 2012). Apparently, this feature 

has been overlooked in previous electron-density studies on 
proteins (Jelsch et al, 2000; Schmidt et al, 2003, 2011; Guillot 
et al, 2008). 

Secondly, we provide a characterization of chemical 
bonding near the active site of HEWL through topological 
analysis of the static electron density. Following two earher 
reports on amino acids (Mondal et al, 2012; Prathapa et al, 
2013), we present here the dynamic electron density of the 
lAM and multipole model of HEWL, along with an analysis of 
their topological properties. Thus, we demonstrate that we can 
compute dynamic densities of structure models of proteins 
which are free of series-termination effects. 



2. Computational details 

X-ray diffraction data of HEWL without substrate, as 
deposited in the PDB (PDB entry 2vbl; Wang et al, 2007; 
Table 1), were employed for structure refinements of the lAM 
and of a database-based multipole model. Prior to the struc- 
ture refinements, a solvent correction of the structure factors 
according to Phillips (1980) was performed (Jiang & Briinger, 
1994). Structure refinements were performed with the 
computer program MoPro (Guillot et al, 2001; Jelsch et al, 
2005) against reflections of the single-crystal X-ray diffraction 
data, employing the strategy described for human aldose 
reductase by Guillot et al (2008). 
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The structure model of HEWL published by Wang et al. 
(2007) (PDB entry 2vbl) was employed as the starting model 
for the lAM refinement. The overall scale factor, the coordi- 
nates of ordered non-H atoms with a temperature factor B of 
<9 and the anisotropic ADPs of these atoms were refined 
alternatingly, thus following the procedure proposed by 
Guillot et al. (2008). Atoms not included in this subsystem 
were not varied within the structure refinement. A refinement 
of all atoms with B < 8 A^ led to marginally worse agreement 
indices, whereas a subsystem of atoms with > 9 A^ would 
contain heavily distributed atoms that are not suitable for 
refinement (Guillot et al, 2008). 

H atoms were geometrically constrained and fixed at bond 
lengths known from X-ray diffraction (Allen, 1986). The 
ADPs of H atoms were constrained to values of 1.2 (for C— H 
and N-H) or 1.5 (for O-H, -CH3 and -NH3) times the 
values of their neighbouring atoms (Guillot et al, 2001; Jelsch 
et al, 2005; Miiller et al, 2006). 

As the final stage of their refinement, Wang et al (2007) 
reported an optimization of the weighting scheme. We have 
not succeeded in obtaining details of this weighting scheme. 
Therefore, we omitted the 293 worst-fitting reflections with 
l^bs — ^caic|/cr(^bs) > 45. A refinement without the employ- 
ment of this criterion resulted in /?f = 8.60% for all reflections. 
The choice of the criterion {I^bs — 4aicl''f(^bs) > 45 reflects the 
possibility that a few reflections might have obtained inaccu- 
rate intensities in the experiment, while by this choice the 
value of is lowered to 8.39% (R^ = 8.54% for all reflec- 
tions). 

Subsequently, C— H, N— H and O— H bonds were elon- 
gated towards the bond-length values known from neutron 
diffraction (Allen & Bruno, 2010). This procedure is necessary 
to counteract the known shortening of bond lengths for 
attached H atoms arising from the shift of the one-electron 
entity, as observable by X-ray diffraction, away from the true 
position of the hydrogen nucleus. More accurate atomic 
positions for H atoms are obtained in this way (Steiner, 2002). 
Refinement of the scale factor resulted in the final lAM. At 
convergence = 8.39% (i?p = 8.54% for aU reflections), 
which is equal to the value obtained with X-ray distances. 

From this model, the static electron density PiAM(r) was 
computed by VMoPro from the MoPro package. This software 
generates the electron density as a superposition of atomic 
electron densities described by the same atomic wavefunctions 
as employed for computation of atomic scattering factors in 
the refinement (Guillot et al, 2001; Jelsch et al, 2005). 

The dynamic lAM density PiAM(r) was constructed by 
superposition of thermally smeared spherical atomic electron 
densities (as defined by the ADPs) by the computer program 
PRIOR (van Smaalen et al, 2003; Mondal et al, 2012). The 
dynamic density of a structure model is basically a Fourier 
map of the structure factors, albeit free from the series- 
termination effects that are necessarily present in bonding 
regions when computing Fourier maps at any experimental 
resolution. ActuaUy, the computer program PRIOR computes 
dynamic electron densities by Fourier inversion of a complete 
set of structure factors up to a resolution of approximately 



dmin = 0 08 A. The structure factors are computed from the 
structure model (coordinates, ADPs, site occupancy factors 
and multipole parameters) employing the same atomic 
wavefunctions as are used in computer programs for multipole 
refinements (Mondal et al, 2012). 

A multipole model was obtained by transfer of multipole 
parameters from the ELMAM2 database (Domagala et al, 
2012) to aU protein atoms of the final JAM. Refinements were 
performed alternatingly of the overall scale factor, the atomic 
coordinates of ordered non-H atoms with a temperature factor 
B < 9 A^ and the ADPs of these atoms, while all multipole 
parameters were kept fixed at their database values. The 
utilization of multipole parameters from the ELMAM2 
database leads to a slightly better fit to the data than the fit 
obtained with the final lAM (Jelsch et al, 1998; Dittrich et al, 
2005; Kalinowski et al, 2007; Johnas et al, 2009; Table 1). The 
resulting multipole model is referred to as the ELMAM2 
model. By superposition of the aspherical atomic electron 
densities based on this model, the static electron density 
PMp(r) was generated by VMoPro in a analogous manner to 
the computation of Piam('')- The computer program PRIOR 
was used to construct the dynamic multipole density Pmp('') by 
superposition of thermally smeared aspherical atomic electron 
densities in an analogous manner as was performed for 

Topological properties of the dynamic electron densities 
were calculated with the computer program EDMA (Palatinus 
et al, 2012), while VMoPro was employed for topological 
analysis of the static densities. Topological properties include 
positions of local maxima of the densities, which are identified 
with positions of atoms, and atomic basins, which are the 
volumes around the maxima that are assigned to each atom. 
The integrated density within a basin provides the atomic 
charge. Furthermore, the density has saddle points between 
pairs of atomic maxima, which are called bond critical points 
(BCPs). According to QTAIM (Bader, 1994), the existence of 
a BCP between two maxima indicates a chemical interaction 
between this pair of atoms. The associated density values 
p(BCP), as well as the values of the second derivatives of the 
density, as provided through eigenvalues of the Hessian matrix 
and Laplacians V^p(BCP), can be used to obtain information 
about the character of the interaction. Kinetic energy densities 
G(BCP) and potential energy densities F(BCP) as well as 
total energy densities //(BCP) and the ratio |F(BCP)|/G(BCP) 
have been calculated according to the approximate formulas 
of Abramov (1997), which allow these quantities to be 
obtained from p(BCP) and VV(BCP). 



3. Results and discussion 

3.1. Global properties of density maps 

All four electron-density maps are positive everywhere, 
with the exception of pM^p(r), which contains a few small 
regions, far away from atoms and bonds, of negative density, 
with a very small minimum value of —5.4 x 10~^ e A~^. Since 
negative electron densities are nonphysical, their presence in 



Acta Cryst (2014). D70, 1136-1146 



Held & van Smaalen • Active site of lysozyme 1139 



research papers 



PM^p(r) indicates a shortcoming of the model, which may be the 
result of using values of the multipole parameters from the 
ELMAM2 database instead of using their true values. Alter- 
natively, negative densities may indicate a non-perfect 
deconvolution of static structure and thermal motion owing to 
the limited resolution of the X-ray diffraction data (Volkov et 
al, 2007). 

While all local maxima of PiAM(r) can be identified with 
atomic positions, PMp(r) contains 38 spurious local maxima in 
low-density regions. VMoPro does not provide this kind of 
information for static electron densities. We again attribute 
their occurrence to incompatibilities between neighbouring 
aspherical atoms from the ELMAM2 database. They do not 
constitute a problem and have not been considered in the 
analysis of bonding features, because they occur far outside 
atoms and bonds (with a distance larger than 3 A). These 
spurious maxima possess a total integrated volume of corre- 
sponding basins of less than 1 % of the volume of the unit cell, 
and contain an integrated number of electrons smaller than 
0.1% of the total number of electrons. 

Within the low-density regions only small differences exist 
between PMp(r) and p^i^^i^), as shown for Glu35 and Asp52 in 
Figs. 2(a), 2(c), 2(e) and 2(g). While the static lAM density 
exhibits nearly spherical contours for p(r) > 3 e A^"*, the static 
multipole density shows pear-shaped distortions of these 
contours pointing into the covalent bonds, thus reflecting the 
description of electron densities within covalent bonds by the 
ELMAM2 model. 



The dynamic electron densities exhibit elliptically distorted 
densities about the atoms, which are especially apparent for 
the O atoms of the carboxyl groups (Figs. 2b, 2d, 2/ and 2h). 
This distinctive feature can be attributed to the anisotropy of 
the atomic displacements. As opposed to small molecules at 
r < 100 K (Mondal et al, 2012), differences between lAM and 
ELMAM2 models are almost obhterated in their dynamic 
electron densities owing to the large values of the ADPs 
(§3.2). 

BCPs of covalent non-hydrogen bonds are found at corre- 
sponding positions in all four electron-density maps studied in 
the present work (Tables 2 and 3). As explained elsewhere, 
thermal motion dramatically diminishes the peak values of 
local maxima, such that some H atoms do not constitute a local 
maximum in dynamic electron densities, and their electron 
density is incorporated into the atomic basin of the X (C, N or 
O) atom to which they are bonded (Hofmann et al, 2007). 
Accordingly, the corresponding Jf— H bonds do not possess 
BCPs (Tables 2 and 3). 

3.2. Atomic displacements and covalent bonding 

A multipole model is a better structure model than the 
lAM. This holds true for a multipole model with multipole 
parameters fixed to values from a database (Afonine et al, 
2004; Dittrich et al, 2005, 2008). While atomic positions might 
differ by very small amounts between the multipole model and 
lAM, this generally means that ADPs in a multipole model are 
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(b) id) (f) (h) 

Figure 2 

Sections of 6 x 6 through the planes of the carboxyl groups of Glu35 and Asp52. (a) pfAM(r) of Glu35; (6) piSviCr) of Glu35; (c) pMp(r) of Glu35; (d) 
p%Hr) of Glu35. (e) pIamCf) of Asp52; (/) p'^SviCr) of Asp52; (g) pMp(r) of Asp52; {h) p%l{i) of Asp52. Contour lines of equal density are given from 0.2 
to 3.0 e A"^ in steps of 0.2 e A"^. 
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Table 2 

Topological properties at bond critical points (BCPs) of Glu35 from static 
and dynamic densities. First line, /o(BCP) (e A~^); second line, V^/o(BCP) 
(eA-^). 
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smaller than ADPs in the corresponding lAM. This feature 
is explained by electron density in covalent bonds being 
mimicked by increased values of ADPs, while it is explicitly 
present in the aspherical scattering factors of the multipole 
model. It is presently observed for HEWL, where the ADPs 
from the lAM are slightly larger than those from the 
ELMAM2 model, although the differences are small. 

The ADPs from the ELMAM2 model have magnitudes 
within the same range of values as possessed by the ADPs 
from the lAM published by Wang et al. (2007). While the 
sample temperature of HEWL was 100 K for the X-ray 
diffraction experiment, the ADPs are considerably larger than 
the ADPs of similar atoms in serine at 100 K and even larger 
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Topological properties at bond critical points (BCPs) of Asp52 from static 
and dynamic densities. First Une, /o(BCP) (e A~^); second line, V^p(BCP) 
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than the ADPs in serine at 298 K (Mondal et al., 2012; 
Fig. 3). 

The discrepancy between expected amplitudes of atomic 
vibrations and values of ADPs can be explained by the 
intrinsic flexibility of proteins (Radivojac et al., 2004; Yuan 
et al, 2005; Weiss, 2007; He et al., 2009). Although triclinic 
HEWL can be regarded as a rather rigid protein with a low 
solvent content and with low B factors compared with other 
HEWL crystal structures, approximately one third of the side 
chains exist in more than one conformation (Wang et al., 
2007). Such a feature can only be partially modelled by 
multiple positions of the atoms, and it will be one cause of the 
large values of the ADPs. Nevertheless, the large values of the 



Acta Cryst. (2014). D70, 1136-1146 



Held & van Smaalen • Active site of lysozyme 1141 



research papers 



ADPs apply to the ordered part of the ELMAM2 model of 
HEWL (Fig. 3). Since at 100 K water and other solvents will 
be frozen, these ADPs will reflect static disorder rather than 
atomic vibrations. 

Whether the observed static disorder represents the flex- 
ibility of HEWL at room temperature or whether its reduced 
magnitude might be caused by the process of cooling cannot 
be determined from the X-ray diffraction data. However, 
ADP values that are larger than the values that are common 
in room-temperature structures of small-molecule crystals 
have the important implication that a meaningful refinement 
of multipole parameters cannot be performed (Dittrich et al. , 
2009; Gatti & Macchi, 2012). This feature is reflected in the R 
values, which are only marginally better for the ELMAM2 
model than for the lAM. This is not the result of the limited 
resolution of the diffraction data or possible poor quality of 
the intensity values (Kalinowski, 2010). As stated above, the 
large ADPs obliterate most if not all of the effects of chemical 
bonding on the electron density and thus on the diffraction. 
Hence, a simultaneous refinement of both ADPs and multi- 
pole parameters cannot be performed in a meaningful way. 

Instead, the refinement of multipole models is only mean- 
ingful with the multipole parameters fixed to values from a 
database (Jelsch et al, 1998; Dittrich et al, 2008; Dominiak et 
al, 2009; Domagala et al, 2012), as performed here for the 
ELMAM2 model of HEWL. Apparently, this feature has been 
overlooked in previous electron-density studies on proteins, in 
which at least part of the multipole parameters have been 
refined (Jelsch et al, 2000; Schmidt et al, 2003, 2011; Guillot et 
al, 2008). 

Despite fixed database values for the multipole parameters, 
the static electron density of the ELMAM2 model does 
provide a good description of covalent bonds. As noticed, the 
BCPs of covalent non-hydrogen bonds are found at 
corresponding positions in PMp(r) and Piam('')- Generally, 
p5li^(BCP) > pSmCBCP) (Tables 2 and 3). Laplacians at BCPs 
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Figure 3 

B factors of Glu35 (red filled squares), Asp52 (red filled circles), B factors 
of ordered glutamic acid (black open squares) and aspartic acid (black 
open circles) residues from HEWL and corresponding B factors for 
serine at 100 K (blue asterisks) and at 298 K (blue diamonds) (Mondal et 
al., 2012). Values are from the particular multipole models. 



of covalent bonds are strongly negative for Pmp('')> while 
values close to zero or even positive values are found for 
VViam(BCP). Large values of p(BCP) and negative Lapla- 
cians of large magnitude are indicative of strong covalent 
bonding according to the QTAIM (Bader, 1994). The differ- 
ences between the two static electron densities at their BCPs 
are the result of chemical bonding, which is solely contained in 
the ELMAM2 model. These features are similar to those 
found for small molecules (Mondal et al , 2012; Prathapa et al , 
2013). They demonstrate the failure of the lAM for proper 
characterization of chemical bonding. 

Primarily, the effect of thermal motion of the atoms on the 
electron density is to lower the local maxima, Le p''''"(r) ^ 
p'""(r) at the local maxima (Mondal et al, 2012). Thermal 
averaging has only a moderate influence on the electron 
densities at BCPs, with p'JX'm(BCP) ~ pSm(BDP) (Tables 2 
and 3; see also Mondal et al, 2012). The fact that p'^p(BCP) ~ 
Piam(BCP) corroborates the observation that large ADPs 
lead to dynamic lAM and multipole densities in which most 
of the effects of chemical bonding are obliterated by the large 
distribution of the atomic positions (Fig. 2). 

The Laplacians at the BCPs of dynamic electron densities 
strongly depend on the values of the ADPs (Mondal et al, 
2012). Here, we find positive values of large magnitude for 
V^PiS'm(BCP) ~ VV^^p(BCP), as opposed to smaU values or 
strongly negative values for VVSm(BCP) and VVmp(BCP). 
These findings are on a par with those for small molecules 
(Mondal et al, 2012). 

3.3. Hydrogen bonds 

Fig. 4 provides a schematic drawing of the hydrogen bonds 
formed by Glu35 and Asp52, as presently identified through 
their BCPs in the static and dynamic electron densities. 

Strynadka & James (1991) reported two hydrogen bonds 
for Glu35. Here, we confirm a hydrogen bond between the 
carboxyl group of Glu35 as an acceptor and the NH group of 
alanine AlallO (0^^_35- ■ -H.llO-N.llO). Instead of forming 
a hydrogen bond to a water molecule, the strongest hydrogen 
bond of Glu35 is between N_35-H_35 and 0_31 of Ala31. 





jl^. /VHd22_59 



Figure 4 

Schematic drawing of the hydrogen bonds formed by Glu35 and Asp52. 
Dashed lines represent hydrogen bonds and green spheres indicate their 
BCPs. 
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Table 4 

Net atomic charges Q of Glu35 and Asp52 from dynamic densities and 
from quantum-chemical calculations (Godjayev et at, 1998). For dynamic 
densities H atoms are included within the atomic basin of their 
neighbouring atom owing to the absence of a local atomic maxima for 
H atoms. 

Atom Qp'^r) Gp&) G quant. 
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N(H) 
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The observed differences between hydrogen bonds in 
different structures of HEWL is probably owing to different 
solvent contents of the crystals used by Strynadka & James 
(1991) and Wang et al. (2007). In addition to these two 
hydrogen bonds, we find six weak hydrogen bonds of the type 
C-H - O (Table 2; Desiraju & Steiner, 2001). Although 
weak, they might be important for defining the stability of 
HEWL and the role of Glu35 in the catalytic process. 

In agreement with Strynadka & James (1991), for Asp52 
we find several hydrogen bonds involving 0*^_52 and 0*^_52, 
although solvent molecules are presently not involved 
(Table 3). Furthermore, we find an important hydrogen bond 
involving N_52— H_52 as a donor. Similar to Glu35, Asp52 is 
involved in six weak C— H- • -O hydrogen bonds, which might 
be important for the stability and role of Asp52 in the catalytic 
process. 

A first indication about the strengths of hydrogen bonds is 
obtained from the values of the densities and Laplacians at the 
BCPs, whereby the static ELMAM2 density appears to be the 
most informative of the four electron densities considered 
here (Tables 2 and 3). The relations between properties at 
BCPs of the four electron densities are in agreement with 
those observed for small molecules (Mondal et al, 2012). 
Values of up to p(BCP) = 0.24 e A^^ indicate that all hydrogen 
bonds are weak or at best of intermediate strength. In 
agreement with results on small molecules, VV(BCP) is small 



and positive, indicating a closed-shell character of the inter- 
action (Bader, 1994). 

An estimate of the character of hydrogen bonds can be 
obtained from the energy densities at the BCPs (§2; Espinosa 
et al, 1998). The ratio |y(BCP)|/G(BCP) < 1 indicates that 
most hydrogen bonds are weak and of mainly electrostatic 
nature (Espinosa et al, 2002). Only 0_31- ■ -N.SS— H_35 
involving Glu35 and 0''^_52- ■ ■H''^^_59-N''^_59 and 
0_44- ■ ■H_52-N_52 involving Asp52 have |y(BCP)|/G(BCP) 
~ 1.1 {i.e. > 1) and total energy densities that are negative 
(Figs. 5 and 6). Accordingly, these three hydrogen bonds have 
a portion of covalency (Cremer & Kraka, 1984) or are of 
mixed covalent/ionic character, but they still are relatively 
weak hydrogen bonds compared with the many hydrogen 
bonds in small molecules. 

3.4. Reaction mechanism 

The atomic charges of the atoms of Glu35 and Asp52 are 
of comparable magnitudes for PiXm(i") and PMp(r) (Table 4), 
although the latter should be considered to be more accurate. 
Integrated atomic charges could not be obtained by VMoPro 
for the static electron densities. As discussed in §3.1, H atoms 
are included in the atomic basins of the non-H atoms to which 
they are covalently bonded. Accordingly, charges are reported 
for NH, CH and CH2 groups for both the experimental and the 
theoretical static electron densities, as indicated in Table 4. 

The atomic charges from PMsit) [and PiAM(r) as well] 
follow chemical expectations, with positively charged C atoms 
and negatively charged O atoms and NH groups. These 
properties are essentially different from the theoretical elec- 
tron densities of Godjayev et al (1998), who proposed a 
negatively charged O' atom of the carboxyl group of Asp52 
along with positively charged O atoms (Table 4). Although not 
impossible in principle, the present experimental electron 
densities lead to the scenario of O atoms of the carboxyl group 
that possess negative charges. This is of importance for the 
proposed catalytical mechanisms of HEWL. According to the 
Phillips mechanism (Phillips, 1966), negative charge of Asp52 
is essential for the enzymatic reaction since it is supposed to 
stabilize the carbenium ion formed after the cleavage of the 
glycosidic bond. 

However, the carboxyl group of Glu35 is not protonated in 
the present crystal structure, which may be owing to pH values 
during crystal growth that do not correspond to those of the 
active enzyme, or it may be the true state of Glu35 in HEWL. 

Since this proton is an essential proton for the Philhps 
mechanism (Phillips, 1966), the chemical characters of Glu35 
and Asp52 of the present structure of HEWL are consistent 
with both the Phillips mechanism (Phillips, 1966) and the 
Koshland mechanism (Koshland, 1953), prohibiting a clear 
approval of one or the other reaction mechanism. 

4. Conclusions 

Wang et al (2007) have noticed that HEWL is rather rigid for 
a crystalhzed protein molecule, but that it nevertheless 
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Figure 5 

Energy densities at BCPs of liydrogen bonds of Glu35. (a) Kinetic energy density G(BCP), (b) potential energy density y(BCP), (c) total energy density 
H{BCV) = G(BCP) + y(BCP) and (d) the ratio |F(BCP)|/G(BCP). Black squares, values from p^|jp(r); red circles, values from pi^AM(r); green triangles, 
values from pMp(r); blue diamonds, values from /ofAM(r). 



contains disordered parts for about 35% of the structure. 
Here, we have found that the ADPs of atoms within ordered 
parts of the structure have magnitudes that are larger than the 
magnitudes of the ADPs of small molecules in their crystal- 
lized state at room temperature, despite the nominal 
temperature of 100 K of HEWL (Fig. 3). Therefore, the ADPs 
of HEWL at 100 K reflect frozen static disorder rather than 
thermal vibrations (§3.2). 

It its generally accepted that multipole refinements provide 
meaningful values for multipole parameters only if the 
magnitudes of the ADPs reflect thermal vibrations of the 
atoms at temperatures of 100 K and below (Gatti & Macchi, 
2012). An important consequence of the observed large ADPs 
is thus to prevent meaningful refinement of the multipole 
model, despite reports of such refinements for some proteins 
in the literature (Jelsch et al, 2000; Schmidt et al., 2003, 2011; 
Guillot et al., 2008). 

Instead, useful information about the redistribution of 
electron density owing to chemical bonding can be obtained 



from refinements of atomic coordinates and ADPs of structure 
models incorporating multipole parameters fixed to values 
from a database of transferable multipole parameters, as is 
presented here for the ELMAM2 model of HEWL (Jelsch et 
al, 1998; Dittrich et al, 2008; Dominiak et al, 2009; Domagala 
et al, 2012). 

We present and analyse static and dynamic electron densi- 
ties based on a multipole model of HEWL. The large ADPs 
are reflected in the dynamic electron densities, such that the 
dynamic lAM and ELMAM2 densities are similar within their 
low-density regions, despite clear differences between the 
corresponding static densities (Tables 2 and 3). The flexibility 
of a protein is visible as distinct differences between static and 
dynamic densities. However, the quantitative interpretation 
of topological properties of dynamic densities is beyond the 
present state of the art of charge-density analysis (Mondal et 
al, 2012). 

Consideration of the static ELMAM2 density has shown 
that electron densities in covalent bonds have similar prop- 



1 1 44 Held & van Smaalen • Active site of lysozyme 



Acta Cryst. (2014). D70, 1136-1146 



research papers 



—I — I — I — I I 



Figure 6 



O 
o 



■ 


Dyn. MP 


o 


Dyn. lAM 




Stat. MP 


O 


Stat. lAM 



A • 

ft 



5 : 



O-40 



-80 



▲ 

O 



A 

o 



8 



8 



o 



■ 


Dyn. MP 


o 


Dyn. lAM 


A 


Stat. MP 


o 


Stat. lAM 



Hydrogen bonds 

(«) 



Hydrogen bonds 

ib) 



—I r- 



—I r- 



—I r- 



8 o 



■ 


Dyn. MP 


o 


Dyn. lAM 


A 


Stat. MP 


o 


Stat. lAM 




Hydrogen bonds 



Hydrogen bonds 

id) 



Energy densities at BCPs of hydrogen bonds of Asp52. (a) Kinetic energy density G(BCP), (b) potential energy density y(BCP), (c) total energy density 
H{BCV) = G(BCP) + y(BCP) and (d) the ratio |F(BCP)|/G(BCP). Black squares, values from p^|jp(r); red circles, values from p'I^amM; green triangles, 
values from pMp(r); blue diamonds, values from PiAM(r)- 



erties to those in small biological molecules. The analysis of 
intramolecular interactions involving the Glu35 and Asp52 
residues has revealed the presence of N— H- ■ -O hydrogen 
bonds of intermediate strength, in agreement with the litera- 
ture. The presence of BCPs in the ELMAM2 density as indi- 
cator for chemical interactions between atoms has revealed 
many weak C— H- ■ O hydrogen bonds, which might be 
important for the stability and function of this protein. 

The topological analysis of the ELMAM2 density has lead 
to negative charges on the carboxylic O atoms of Asp52, thus 
supporting the Phillips mechanism for the catalytic activity of 
HEWL (Phillips, 1966). However, the deprotonated state of 
the carboxyl group of Glu35 makes the present electron 
densities in agreement with both the Koshland and Phillips 
mechanisms (Koshland, 1953; Phillips, 1966), and they cannot 
be distinguished here. Electron densities of one or more 
complexes of HEWL with a substrate or inhibitor and an 
electron density of HEWL with Glu35 in the protonated state 
might allow a unique assignment of the mechanism and might 



lead to a more detailed characterization of the catalytic 
mechanism. 

Financial support of the present work by the German 
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